Entangled Webs: Tight Bound for Symmetric Sharing of Entanglement 
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Quantum entanglement cannot be unlimitedly shared among arbitrary number of qubits. Larger 
the number of entangled pairs in an iV-qubit system, smaller the degree of bi-partite entanglement 
is. We analyze a system of N qubits in which an arbitrary pair of particles is entangled. We show 
that the maximum degree of entanglement (measured in the concurrence) between any pair of qubits 
is 2/N. This tight bound can be achieved when the qubits are prepared in a pure symmetric (with 
respect to permutations) state with just one qubit in the basis state |0) and the others in the basis 
state |1). 



Schrodinger fjj] has identified quantum entanglement 
as the key ingredient in the paradigm of quantum me- 
chanics. Throughout the whole history of modern quan- 
tum mechanics the mystery of quantum entanglement 
puzzled generations of physicists g . On the other hand 
in the last decade the entanglement has been recognized 
as an important resource for quantum information pro- 
cessing. In particular, quantum computation J 3|, quan- 
tum teleportation quantum dense coding certain 
types of quantum key distributions || and quantum se- 
cret sharing protocols 0, are rooted in the existence of 
quantum entanglement. 

In spite of all the progress in the understanding of 
the nature of quantum entanglement there are still open 
questions which have to be answered. In particular, it 
is not clear yet how to quantify uniquely the degree of 
entanglement [p|-|T2|] , or how to specify the inseparability 
conditions for bi-partite multi-level systems (qudits) [ p"3|| . 
Further problem which waits for a thorough illumination 
is the multiparticle entanglement fll4j . There are sev- 
eral aspects of quantum multiparticle correlations. For 
instance, it is the investigation of intrinsic n-party entan- 
glement (i.e, generalizations of the GHZ state p5| ). An- 
other aspect of the multiparticle entanglement is that in 
contrast to classical correlation it cannot freely be shared 
among many objects. In particular, Coffman et al. [lfj] 
have studied a set of three qubits A, B and C. It has 
been shown that the sum of the entanglement (measured 
in terms of the tangle ) between the particles AB and 
the particles AC is smaller or equal to the entanglement 
between the particle A and the subsystem BC. Woot- 
tcrs has considered an infinite collection of qubits 
arranged in an open line, such that every pair of nearest 
neighbors is entangled. In this translationally invariant 
entangled chain the maximal closest-neighbor (bi-partite) 
entanglement (measured in the concurrence) is bounded 
by the value l/y2 (it is not known whether this bound 
is achievable) fL7|) . 

In this Letter, we consider a finite system of N qubits 
in which each pair out of N(N — l)/2 possible pairs is 
entangled. We show that the maximal possible bi-partite 
concurrence in this case is equal to 2/N . The derivation 
of this tight bound on the concurrence is the main result 



of our Letter. 

The problem is formally posed as follows. Suppose that 
N qubits, indexed by I = 1, 2, . . . , N, are in the state pat, 
and choose a basis {|1), |0)} for each qubit. We require 
that the matrix form of the marginal density operator 
pu> for a pair of qubits I and V , represented in the cho- 
sen basis, is independent of the choice of I and V . Note 
that this requirement is satisfied if p^ is invariant un- 
der any permutation of qubits. The question is to find 
the maximum degree of entanglement between a pair of 
qubits. 

It is convenient to suppose that each qubit is a spin- 
1/2 particle with the spin operator s^ l '(l = 1,...N). 
The Hilbert space of the subsystem composed of qubits 
1 and 2 is a direct sum of the subspaces for the 
total spin and 1, with the projectors Pq and P\ 
onto each subspace, respectively. Under the condition 
Pu = P2i(= p), we have PopPi — PipPo = 0, since 
these operators change their signs under the permuta- 
tion of the two qubits. Let us define irreducible tensors 
Tjg of rank k = 0,1,2 and components q such that 



T$% = E m w<A,9|l.m;l s mO(ELi^)(Ej=iC ; ), 
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where = ± %W)/y/2, 4 

(k, q\l, to; 1, to') is the Clebsh-Gordon coefficient for 
forming a total spin k state from two spin-1 particles. 
The spin-1 part of the density operator p can be expanded 

by as P\pP\ — J2k q a M-^2,g ' an d t ne coefficients 

ak,q are obtained by the relation ak,qTr(T 2 _ T 2 k q) = 



Ti{f 2 (k l q p) = (T 2 {K f q ), where we denote Tr(. . . p N ) as 
(...). From the symmetry of p^, we have (Tj^ 



=.(*) 



, Nql - (N(N - 1)/2)(T 2;9 
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Pn given, it is convenient to choose the x,y, and z axis 
as the principal axes for the tensor of the second-order 
correlation for the total spin S = YaLt. ^ °f t ne N 
qubits, namely, (S, l S v + S^S^) = 2^5^ with S* = (S^) 
(here and henceforth, suffices p and v represent x,y,z). 
Matrix elements for p then takes a simple form on the 

basis {|TT) - |||>, |TT> + III), III) + |U>, |U) - |U» as 
follows, 
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where we have introduced nonnegative parameters Aq = 

N(N ^N-ij S ) and A v = ^(Af-i)' - Aq, that satisfy 
A x + Ay + A z + A = N. 

As a measure of entanglement between the two qubits, 
we use the 'concurrence' introduced by Hill and Woot- 
ters Q. The concurrence C can be calculated as fol- 
lows. Let p be the time reversal of p, which is obtained 
by changing the sign of spin (S) in the expression (Q). 
The eigenvalues of pp are all real and non-negative, and 
let the square roots of those be h, h, h, and U in 
decreasing order. The concurrence C is then given by 
C = maxjZi — Z2 — Z3 — Lt, 0}. In the present case, one of 
the eigenvalues of pp is (A /N) 2 . Let us denote the other 
three as (Ai/iV) 2 , (A 2 /iV) 2 , and (A 3 /iV) 2 in decreasing 
order, and introduce parameters (3 = X\ + A 2 + A3 and 
7 = Ai — A2 — A3. The concurrence is then given by 
C = max{(7 — A Q )/N, (A — [3)/N,Q}, where allowances 
are made for the possible order of Aq and \\. 

In the following, we first fix the parameters S 2 (hence 
A^ and Aq), an d maximize 7 with respect to (X, Y, Z) = 
((Sz) 2 , (S y ) 2 , (S z ) 2 ). We then move S 2 to obtain the 
global maximum of (7 — Aq)/N, that turns out to be, 
as we shall see, the maximum of the concurrence. For 
simplicity, we assume that S 2 > S 2 > S 2 (hence A z < 
A y < A x ). The states with some parameters equal will 
be considered as the limiting cases. 

There are two simple bounds for the allowed values of 
(X, Y, Z) . One is a necessary condition for p to be phys- 
ical. The eigenvalues for p must be nonnegative, and the 
boundary is given by the surface that satisfies det(p) = 0. 
Calculating from (|l|), this surface turns out to be a plane, 
and the condition for (X, Y, Z) is 



j a = A x A y A z - A X X - A y Y - A Z Z > 0. 



(2) 



The other one is a requirement necessary for the 
spin correlations. From the inequality (Eji^i)^ — 
(S l ))/S 2 } 2 ) > 0, we obtain 



fs = l 
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S 7 



Y 



Z 



s 2 



s 2 



> 0. 



(3) 



Any physical state thus falls in the region V, that is de- 
fined by f A > 0, f s > 0, and X,Y,Z>0 

The relations that connect (X, Y, Z) and Ai are ob- 
tained by expanding det[pp ~ [\ 2 /N 2 )I] = Y[-(\ 2 - A 2 ) 
and equating the coefficients of A™. There are three in- 
dependent equations, and it is convenient to take the 
following set: 

fo = A\ + A 2 y + A\ - 2X - 2Y - 2Z 
= \\ + \\ + A 2 



(4) 



4Y AZ\ 

Ib = —B X B V B Z ( B x + By + B z — —J 



B, 



B,, 



— Ai + A 2 + A 3 — 2{\ 1 X 2 + A 2 A 3 + AgAj), 

where B^ = A x +A y +A z ~2A fl . The relation A x > A y > 
A z > implies that B z > B y > \B X \. Note that the sign 
of Jb is the same as that of 7, since we can factorize as 
f B = + A 3 - A 2 )(Ai + A 2 - A 3 ). 

Let W be the region defined by Ja > and /s > 0. 
As a function of {X, Y, Z), 7 is continuous in the re- 
gion W including the boundaries. The gradient V7 = 
(dj/dX, dj/dY, dj/dZ)ca.n formally be obtained by us- 
ing the three relations ([|). The result is 



c?7 
dX 



( 1 + B y )( 7 + B z ) 



(5) 



where n = 2(Ai — A 2 )(Ai — As)(A 2 + A3) > 0, and the 
other two are obtained by the cyclic exchange. At the 
inner points of W, Ai > A 2 since 7 > 0, and A 3 > 
since /a > 0. The parameter k is hence positive, and 
the gradient V7 exists. Since B z > B y > 0, we have 
dj/dX > for the inner points of W. 

The geometry of W and V is derived as follows. Let 
the points on which the planes /b = 0, Ja = 0, and 
fs = intersect X-axis be Pbx, Pax, and Psx, re- 
spectively, and denote other six points on Y- and Z-axis 
similarly. The relation A x A y - (B x + B y + B Z )B Z /A = 
B x B y /A, and the ones obtained by the cyclic exchange 
of {x,y,z}, tells us the following. When B x > 0, the 
triangle ttb = PbxPbyPbz does not intersect with 
ir a = Pax Pay Paz, and lies closer to the origin. W is 
the sandwiched region of the two triangles. The triangle 
■Kg = PsxPsyPsz may intersect with tta and ttb or not. 
When B x < 0, fs > (hence 7 > 0) is satisfied by all the 
points that satisfy Ja > 0. When B x — 0, 7 = for the 
points on yZ-plane, and 7 > for X > 0. Combining 
these observation with dj/dX > 0, we conclude that in 
V, 7 takes its maximum on points on the boundaries tta 
or ns, and never on the inner points. 

To determine the behavior of 7 on the boundaries, we 
first derive the value of 7 on the axes explicitly. For the 
points on Z-axis and satisfy (^|), the roots {Ai} are found 



to be {A 



(A x + A v ) 2 - A{S Z ) 2 ± (A x - Ay) 
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The expression for 7 depends on which is the largest root. 
The roots for the points on the other two axes are simi- 
larly obtained. Applying these for the vertices of tta, we 
have j{Pax) - B y and j(P A y) = i{Paz) = \B X \. On 
the vertices of 7rs, 7 is well-defined only when they are 
in the region W . When P$z is in W, 



< i(A x A y - S 2 ) 
= (A Z + A Q )(A Z - 



A - 2) - [A x ~ A y ) 2 



(6) 



and hence A z > 2 — Aq . Then 7 is written as 
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l(Psz) = A Z - y/ (A x + A y y - 4S 2 

= A a - V(^z + ^o)(A 2 +A -2), (7) 

This is a decreasing function of A z and noting S* 2 > 0, 
we have 2 - A > 7(Psz) > N/(N -1)-A = j m . Simi- 
larly, j(Psx) > 7m and 7(Psy) > 7m if those points are 
in W. 

Next, we will evaluate the gradient on the bound- 
aries. On the boundary tta, k is not always positive, 
but k = is possible only for the following two cases. 

(a) Ai = A 2 . This occurs only when /a = /b = 0, 
or equivalently, on the segment PayPaz when B x = 0. 

(b) A2 = A 3 = 0. This occurs if and only if _/a = 
and (/o) 2 — /b =0. The intersection of /a = and 
(/o) 2 — /b = is a parabola. The point Pq(A 2 (A x — 
A y )/{A X - A z ),0,Al(A v - A Z )/(A X - A z )) is on this 
parabola. At Pq, u = (A x , A y , A z ) is normal to Ja = 0, 
and v = (—A x By, A y B y — 2A x A Zl —A z B y ) is normal to 
(/o) 2 — Jb = 0. Since mxd has a vanishing Y component, 
the parabola is tangent to JTZ-plane at Pq. Another 
point {A 2 y (A x - A Z )/{A X - A y ),-A 2 x {A y - A Z )/{A X - 
A y ),0), on which Y < 0, is also on the parabola. We 
thus conclude that the gradient exists on tta except for 
the segment PayPaz and Pq. 

Let us define three particular directions on tta as 
qv x = (A y A z ,-A z A x ,0), q zx = (A y A z ,0,-A x A y ), and 
q yz = (0, — A Z A X , A x A y ). The differential coefficients of 
7 for these directions are calculated to be q yx ■ V7 = 
A Z (B 2 Z - ~f 2 ){A x - Ay)/ K , q zx ■ V 7 = A y (B 2 y - 7 2 )(A C - 
A z )/k, and qv> • V 7 = A,( 7 2 - B 2 x ){A y - A z )/k. Not- 
ing that j(Pax) = B y and j(Pay) = i{Paz) = \B x \, we 
conclude that q yx ■ V7 > 0, q zx ■ V7 > and q yz ■ V7 > 
on tta except for on the segment PayPaz and PqPax- 
On PqPax, Q zx • V7 = 0, so that 7 is constant. 

Similarly, define directions on its as p xy = (-S 2 , S 2 ,0) 
p yz = (0,-S 2 ,S 2 z ) and p xz = (S 2 x ,0 7 S 2 ). Let ir sv 
be the intersection of 7r5 and W. On nsVi we have 
pfv • V 7 = (7 + S,)( 7 - lm ){S 2 y - S 2 x )/n , p xz ■ V 7 = 
{ 1 + B v ){ 1 - lm ){Sl-S 2 x )/n, andp^-V 7 = (7+Bx)(7- 
lm){S 2 z - S?)/k- When < 0, 7 (P A y) = l{P A z) = 
\B X \ and (^) implies that 7 = |B X | on yZ-plane, and 
7 + B x > for X > 0. Since we have seen that 7 > j m 
on the vertices on irsv, P xz ■ V7 > and p 311 ' • V7 > 



When B a > 2 - Ao, we have tp < 4 - A (2 - A ) < 4. 
We thus conclude that 7 < 2 for all values of B y . When 
N < 6, 7 approaches 2 only in the limit of Aq — > 0, 
_B a — > 2 , — > 2, and A x — > 2. Pi approaches Pax 
in this limit, so the limit is the same as in the case ii). 
When N > 6, 7 -> 2 only in the limit of A -> 0, P y -> 2, 



everywhere on ttsv, and p yz • V7 > on irsv except for 
segment PsyPsz- 

Now we are in a position to find the maximum of 7. 
We must consider the following four cases separately (see 
Fig. [J). 

i) A x A y > S 2 Z and A y A z > S 2 . In this case, A z > 
2 — Aq is necessary. 7 takes its maximum on P$z and 
the value is given by (Q) and 7 < 2 — A . 7 approaches 
2 only in the limit of Aq — > 0, A z — * 2, Pj; — > 2, and 
B y — > 2. This limit can be taken only if iV > 6, since 
Ay > A z must hold in the limit. When N > 6, Psz is 
the only point that attains 7 = 2 since p yz ■ V7 > and 
p xz • V7 > still hold in the limit. When N = 6, 7 = 2 
everywhere on tta , but these states are equivalent in the 
sense that they are related by the rotation of the whole 
system. 

ii) A x A y < S 2 and A y A z < S 2 . In this case, from 
the relation < 4(5^ - A y A z ) = ~(B y + A - 2){B y + 
A ) - (B z - By)(B y + A - 1), By < 2 - A is neces- 
sary. 7 takes its maximum on PqPax and the value is 
7 = B y < 2 — Aq. 7 approaches 2 only in the limit of 
Aq — ► 0, A x — > 2, P y — > 2, and P 2 — > 2. This limit can 
be taken only if JV" < 6, since A x > A y must hold in the 
limit. When TV < 6, Pq coincides with Pax in the limit 
and Pax is thus the only point that attains 7 = 2. When 
TV = 6, the limit is the same as in the case i). 

iii) A x A y > S 2 and A y A z < S 2 . The maximum is the 
larger of j(Pax) and j(Psz)- Depending on N, one of 
them or both can approach 2. The limit is the same as 
described in the cases i) and ii). 

iv) A x A y < S 2 and A y A z > S 2 . In this case, from the 
relation C(Pax) > C(Psx), we have B y > j m . We 
also have Aq < 2 since (2 - A )(B y + A Q + 2A Z ) = 
4(5 2 - A x A y ) + 2A Z + B y (N -2 + 4S 2 )/(N - 1) > 0. 
7 takes its maximum at Pi, that is the intersection of 
Paz Pax and P S zPsx- Since A 3 = at Pi, 7 = A x - A 2 
and (3 = \\ + A 2 . Then we have 1 + 7 2 = 2/ and 
/? 2 7 2 = /s. This implies that 7 2 is the smaller of the two 
roots t = t a , tp of the equation t 2 — 2fgt + /b = 0. The 
coordinates of Pi can explicitly be obtained by solving 
Ja = and fs = with Y" = 0. Substituting the re- 
sult into /o and fg in the equation of t, we finally obtain 
i Q = P 2 and 



(8) 



B x — > 2, A z — > 2, and Pi — > Pgz. This is the same limit 
as in the case i). When N — 6, 7 — > 2 in the limit of 
Ao — > and all A M approaching 2. 

Finally, we show that Aq — (3 is smaller than 2. In 
PF, Jq takes its minimum on Paz, and its value is 
(A x - A y ) 2 + A 2 . We thus have A - j3 < A - ^% < 



4A X A Z 

(By - lm )(N - 1) 

(2-Aq) 2 + 2A B " 



Aq(A -2) 

-Aq-2 2(B y + A - 2){B y + N) + (B z - B y )(B y 



B^r 



B u 



(N-l)(B y - lm ) 
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2S 2 - 2S*)/(N - 1) < N/(N - 1) < 2 



A - A z = (N ■ 
for N > 3. 

Combining all cases, we reach a conclusion that the 
maximum value of the concurrence is 2/N, and this value 
is reached only by the state satisfying the constraints 
(S*> = (N/2 - l) 2 , (S z > = AT/2 - 1, (S 2 ) = (§*) = 
(3N - 2)/4, and 0a,) = (S'y) = 0, if z-axis is suitably 
chosen. Such a state exists - it is the eigenstate of S z 
with eigenvalue N/2 — 1 with total spin N/2. This state is 
equally- weighted in-phase superposition of any one qubit 
being in the state |0) and the other N — 1 qubits in the 
state This is a permutationally invariant pure state 
that highly (in the order of N) breaks the symmetry be- 
tween the basis states |0) and |1). 

It is worth noting that the optimal state for our prob- 
lem was found to be a pure state. This is not trivial 
because the symmetry is required for density operators. 
Obviously, possible mixed states span a larger Hilbcrt 
space than that spanned by possible pure states. 

For entangled chains, the best state so far is a mixed 
state |0. It will be interesting to investigate entangled 
chains composed of a finite number of qubits. Specifi- 
cally, instead of an infinite line of qubits one can consider 
an entangled loop composed of even number of qubits 
2N (here we assume N > 1 because for TV = 1 we 
would have a two-qubit "loop" which is maximally en- 
tangled when prepared in a Bell state, in this case the 
concurrence is equal to unity). One can find a pure 2N- 
qubit rotationally-invariant state for which the closest- 
neighbor bi-partite entanglement has the concurrence 
C = (2 + 2 N - 2 )/(2 + 2 N ). In the limit N -► oo this 
concurrence is Coo = 1/4- We conjecture that this is the 
maximal value of the concurrence for pure infinite entan- 
gled chains. This means that impure states in the limit of 
infinite number of particles may attain higher bi-partite 
entanglement than pure states. 
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FIG. 1. Allowed regions for the parameters {X, Y, Z}. The 
shaded surfaces are boundaries for the physically allowed 
states. The arrows show the directions in which the con- 



currence increases. 
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